Preprint typeset in JHEP style - HYPER VERSION 



DAMTP-2007-66 



Self-gravitating radiation in AdSd 



Vladislav Vaganov 

D.A.M.T.P., Centre for Mathematical Sciences, University oj Cambridge, Wilberforce 
Road, Cambridge CB3 OWA, U.K. 



E-mail: vv205(§cam . ac . uk 



Abstract: We study spherically symmetric equilibrium configurations of self-gravitating 
massless thermal radiation in asymptotically anti-de Sitter space. In d = 4, it was shown 
by Page and Phillips that there is a maximum red-shifted temperature, maximum mass 
and maximum entropy. For higher central densities, the temperature, mass and entropy 
undergo an infinite series of damped oscillations, corresponding to unstable configurations. 
We extend this work to all dimensions d >2>. We find that in 4 < d < 10, the behaviour 
is similar to the d = 4 case. In d > 11, the temperature, mass and entropy are monotonic 
functions of the central density, asymptoting to their maxima as the central density goes 
to infinity. In d = 3, an exact solution is given by a slice of the AdS C-metric. 
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1. Introduction 

It is well known that the canonical ensemble is not defined in asymptotically flat space. 
This is because having thermal radiation at constant temperature at infinity is not com- 
patible with asymptotic flatness. One can avoid this problem by enclosing the system 
in a box, which is unphysical Q, or by working in asymptotically anti-de Sitter (AdS) 
space. In AdS the canonical ensemble is given by a Euclidean path integral over all matter 
fields and metrics which tend asymptotically respectively to zero and to AdS identified 
with period (3 = T^^ in imaginary time, where T is the red-shifted temperature Q. In 
their seminal paper. Hawking and Page studied this path integral in the semiclassical 
approximation, where it is dominated by classical solutions to the Einstein equations with 
these boundary conditions. Their work was later generalized to d dimensions Q. For 
T < Tq = \/ {d — l){d — 3)/2ttI, where / ^ Ipianck is the AdS length, there are no black 
hole solutions; thermal radiation, that is, AdS periodically identified in imaginary time, 
is the only admissible phase. For Tq < T < Thp = {d — 2)/27r/, thermal radiation is the 
preferred phase, although black holes may form and evaporate from time to time as a result 
of fluctuations. For T > Thp, the configuration with a large AdS black hole in equilibrium 
with thermal radiation has a lower free energy than thermal radiation alone and it there- 
fore dominates the path integral. The point T = Thp marks a first order phase transition, 
known as the Hawking-Page transition, between two topologically distinct manifolds. The 
small AdS black hole has negative specific heat and is never a dominant phase, but it serves 
as a bounce mediating the tunneling amplitude between thermal radiation and the large 
AdS black hole §. 
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Thus far the gravitational effect of thermal radiation has been neglected, and therefore 
this picture will undergo quantum corrections when the back-reaction is taken into account. 
These correspond to the one and higher loop terms in the path integral. An exact solution 
describing a black hole in equilibrium with thermal radiation remains elusive, although 
this problem has recently been conjectured to be equivalent to the quest for a localized 
black hole on a brane |^] , which is a purely classical problem in GR. For thermal radiation 
alone, however, a useful model may be obtained by treating the radiation g perfect 
fluid with equation of state p = p/{d — 1), where the energy density is related to the local 
temperature by the Stefan-Boltzmann law, p cx Tf^^. This amounts to taking just the 



leading term in the high-temperature expansion of the one- loop effective action |10| and 
is a good approximation if the red-shifted temperature is much greater than the inverse 
of a characteristic curvature radius of the spacetime. For a given central density, one may 
solve the Einstein equations with a A term for the static spherically symmetric equilibrium 
configurations of self-gravitating thermal radiation. In d = 4, this calculation was carried 
out by Page and Phillips Q . Their key finding was that there exist locally stable radiation 
configurations all the way up to a maximum red-shifted temperature T^ax ^ Thp, above 
which there are no solutions.^ Thus for temperatures T > Tmax, the only possible phase is 
a large AdS black hole in equilibrium with thermal radiation . There is also a maximum 
mass and maximum entropy configuration occurring at a higher central density than the 
maximum temperature configuration Q. Beyond their peaks, the temperature, mass and 
entropy undergo an infinite series of damped oscillations; configurations in this range are 
unstable, although the precise onset of instability depends on the ensemble The purpose 
of this paper is to generalize the analysis of to d dimensions. 

The organization is as follows. In section |2| we write down the governing equations 
for this system, which are a special case of the Tolman-Oppenheimer-Volkoff equations. 
The solution space is two-dimensional, but if we impose regular boundary conditions at 
the origin we get a one-parameter set of equilibrium configurations labeled by the central 
density. These solutions must be found numerically, although useful asymptotic expressions 
may be obtained when the radius is much smaller or much larger than the AdS length. We 
show how given a solution, one may compute its red-shifted temperature, mass and entropy. 
We have calculated these for a range of central densities in each dimension. The results 
are presented in section ^, where we see that in each dimension d > 3 there is a maximum 
red-shifted temperature, mass and entropy. However, whereas in 4 < d < 10 the behaviour 
is similar to the d = 4 case, in d > 11 the behaviour is different: the temperature, mass and 
entropy increase monotonically with the central density, with no oscillations. In section ^ 
we briefly address the range of validity of the perfect fluid approximation. We conclude 
in section |5| with a discussion of the broader relevance of this work. Details of the d = 3 
case are contained in Appendix ^ Appendix ^ contains an analysis of the scale-invariant 
A = equations. 



""^This temperature had previously been estimated in M: there it is known as T2. 
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2. Perfect fluid radiation in asymptotically anti-de Sitter space 



We want to solve the Einstein equations with radiation gas perfect fluid source for a static 
spherically symmetric metric in d spacetime dimensions, with asymptotically AdS bound- 
ary conditions. We will work in the system of units in which c = h = ks = i and k = SttGn. 
The following is mostly based on Q extended to d dimensions. The metric may be written 
as 

ds^ = -e'^^Vdf + V'^dr^ + r^dnj_2 
/ s / 2K,m(r) \ 

where dQ'^_2 is the metric on the unit (d — 2)-sphere, which has volume E. The AdS length 

1 /2 

/ = {—{d — l){d — 2)/2A) ' , where A < is the cosmological constant. Written in this 
way, 

lim m(r) = M (2.2) 

r— >oo 

is the total mass of the configuration, and the time coordinate t may be normalized such 
that 

lim tlj{r) =0. (2.3) 

r— >oo 

The energy-momentum tensor is taken to be that of a perfect fluid with equation of state 
p = p/{d — 1), where p is the isotropic pressure. When the fluid is at rest with respect to 
the static frame defined by the Killing vector d/dt, this takes the form 

where the rest frame energy density is related to the local Tolman temperature Tioc through 
the Stefan-Boltzmann law, 

p = aTt . (2.5) 

The dimension-dependent constant a is proportional to the effective number of degrees of 
freedom (see section ^ . 

The Einstein equations G^y + I^g^u = kT^u for the metric ( ^.1| ) and the energy mo- 
mentum tensor (|2.4D yield the following system of ODE's: 

V d-2 

dr 

dp _ pd[{d-l){d-3)Km/J: + Kr'^-'^p+{d-l){d~2y-^/P] 

d^ ~ {d - l)id - 2)rd-^V ' ^^'^^ 

which are a case of the Tolman-Oppenheimer-Volkoff equations in d dimensions. There is a 
two-parameter family of solutions, but imposing regular boundary conditions at the origin, 

piO) = Pc m(0) = 0, (2.7) 

yields a one-parameter set of regular equilibrium configurations m{r), p{r) labeled by the 
central density pc > 0. In general ( |2.6| ) must be integrated numerically. Because these 
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equations are singular at r = 0, the solution must be cut off at r = e <C /, where we set 
p(e) = pc and m(e) = T,pce'^~^ /{d — 1). The physical answers are recovered as e is taken to 
zero. In practice, one takes a value of e small enough for the required level of accuracy.^ 

Once m(r) and p{r) (and hence the metric function V) are determined, we turn to the 
contracted Bianchi identity (equivalently, conservation of T^,y), which reads 

-f+pd + =0. (2.8) 



dr \dr 2 dr ^ 

The red-shifted temperature is defined as 

T^Tioc\9u\'^'' = Tioce'^V'/\ (2.9) 

Equation ( p.8[ ) implies that T is constant, as it must be for an equilibrium configuration. 
It follows from condition (|2.3| ) that for large r the solutions behave as p ~ y-'^/^ ^ j.-d,_ 
also follows that given a solution m{r), p{r) the red-shifted temperature may be computed 
by taking the limit 

T = lim 

r— >oo 

Equation (|2.9| ) may then be used to determine ?/'(r). Note that for the regular solutions we 
are considering, equations ( ^.61) and ( |2.8D imply dm/dr > 0, dp/dr < and dip/dr > for 
all r > 0. Pure ^dS* in global coordinates corresponds to pc = and has m = p = ip = Q 
for all r > 0. 

One may compute the entropy S of an equilibrium configuration by integrating the 
local entropy density, 

aTf-^ = -A_a^/dp{d-i)/d (2.11) 



d-l d-l 
over the proper spatial volume: 

S = ^ai/'^S r pid-i)/dy-i/2^d-2 (212) 

If the mass ( |2.2| ) and temperature ( 2.10| ) are known for a range of central densities, the 
entropy may also be obtained by integrating the first law of thermodynamics, dM = TdS, 
starting from the Pc = configuration which has M = T = S = 0. 

We now go on to consider certain limits of the solutions. Observe that equations ( |2.6| ) 
are invariant IH under the scaling 

r^Cr, p^C-^p, m^C^-^m, I ^ CI , (2.13) 

where C G M"*". When A = (/ = cxo) this is an exact scale invariance [^. In this case 
there exists a special self-similar solution^ given by 

p = ar-\ m = aS/-V(t^ - 3) where ^ ^ ^(d - l)(d - 3)^-1 ^ 

d"^ — d + 2 



^It is possible to regularize the Tolman-Oppenheimer-Volkoff equations and recast them as a 3- 
dimensional regular dynamical system on a compact state space This has a number of advantages, 

but is not essential in our analysis, see, however, Appendix ^. 

^By this we mean the self-similarity of the geometry as captured by the existence of a proper homothetic 
vector [^. This does not necessarily follow from the self-similarity of the matter fields as expressed by the 
scaling relation (2.13) with A = 0. Indeed, other solutions of the A = equations do not admit a proper 
homothetic vector. 
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and we take d > 4 (see Appendix ^ for the d = 3 case). This solution is singular at r = 0. 
It is the limit as the central density goes to infinity of the A = solutions with a regular 
origin. The latter may be computed numerically, see, for example, [Q, 17] for d = 4. It 
is best formulated (Appendix |^) using variables invariant under the scale transformation 
(2.13), whence equations ( |2.6| ) may be reduced to a plane autonomous system, from which 



one may derive various qualitative features and asymptotic approximations |Q, 17]. For 
large r, regular solutions will approach the self-similar solution ( |2.14| ), which has infinite 
mass and is thus not asymptotically flat. To obtain finite mass solutions in the A = case, 
one must confine the radiation to an unphysical box. If one does so, the thermodynamics 
(red-shifted temperature, mass and entropy as functions of the central density) of these 
configurations will be qualitatively similar to the asymptotically AdS case studied here. 

For fixed A < 0, equations ( p.6D are no longer scale invariant. However, for r ^ I, the 
terms involving A may be neglected, and the solutions approach the regular A = solution 
with the same value of pc- In the limit as the central density goes to infinity, they approach 
the singular self-similar solution ( p. 14 ) for r <^ I. On the other hand, for large r, p ^ r~'^, 
and the first of equations ( |2.6| ) implies that the mass of these solutions is finite. This is due 
to the confining nature of the gravitational potential in anti-de Sitter space, which acts as 
a box of finite volume. 

3. Results 

For each dimension, we computed a number of solutions for a range of values of pc, and 
the red-shifted temperature, mass and entropy for each solution. A combination of Runge- 



Kutta methods were used to integrate equations (2^) starting from r = e and terminating 



at r = ri ^ I, and then equations (|2.1C| ), (2.2) and ( 2.12 ) with infinity replaced by ri 



were used to evaluate the temperature, mass and entropy. The configuration with infinite 



central density was computed by imposing ( 2.14 ) as a boundary condition at r = e. The 



errors coming from neglecting the region ri < r < oo may be approximated by using 
asymptotic expansions ]p. The accuracy of the calculation depends on the values of e, ri 
and the error tolerance of the routine. A range of different values were used to achieve four- 
significant-figure accuracy, where greater precision was required at higher central densities 
and in higher dimensions. We work in terms of the dimensionless quantities 

Pc = K^Pc 

We found that for all d > 3 there exists a maximum red-shifted temperature, Tmax, a max- 
imum mass Mmax and a maximum entropy Smax, where the maximum entropy configura- 
tion always coincides with the maximum mass configuration. The maxima of temperature, 
mass and entropy, along with the corresponding values of pc, are presented in Table || for 
3 < d < 13. For comparison, we also give the ratio of Mmax to the mass Mhp = (d — 2)S 
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of the large AdS black hole at the Hawking-Page temperature: this is simply to illustrate 
the similar dimensional dependence of these quantities (see section ^ for some physical 
comparisons). Figs. 1-3 show how the dimensionless temperature, mass and entropy vary 
with the central density for three representative cases: d = 4, d = 11 and d = 3. 

















d 


In pc 


f 


In pc 


M 


■^'-^11 iUjJL 1 -^'-^ ri 1 


c 


3 





2/3 ^ 0.6667 


ln4 » 1.386 


7r/3 ^ 1.047 


1/3 ^ 0.3333 


42/37r « 7.916 


4 


1.353 


0.9479 


2.018 


11.56 


0.4598 


15.22 


5 


2.349 


1.081 


2.650 


22.98 


0.3880 


25.95 


6 


3.234 


1.145 


3.354 


37.76 


0.3586 


39.22 


7 


4.119 


1.177 


4.158 


53.39 


0.3444 


52.78 


8 


5.109 


1.190 


5.119 


66.55 


0.3354 


63.81 


9 


6.390 


1.195 


6.392 


74.50 


0.3278 


70.08 


10 


8.590 


1.195 


8.590 


76.08 


0.3204 


70.72 


11 


oo 


1.192 


oo 


71.81 


0.3129 


66.25 


12 


oo 


1.188 


oo 


63.31 


0.3055 


58.14 


13 


oo 


1.183 


oo 


52.57 


0.2983 


48.15 



Table 1: Maximum red-shifted temperature, maximum mass and maximum entropy configurations 
of spherically symmetric thermal radiation in AdS in various dimensions. The quantity Mhp — 
{d — 2)1] is the mass of the large AdS black hole at the Hawking-Page temperature. 
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Figure 1: d = 4: lOT (solid line), M 
(dashed line) and S (dotted line) vs. In pc 
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Figure 2: d = 11: lO^f (solid line), M 
(dashed line) and S (dotted line) vs. In pc 
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Figure 3: d = 3: lOT (solid line), M 
(dashed line) and S (dotted line) vs. Inpc 
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-3.2 -2.4 -1.6 -0.8 0.0 0.8 

Figure 4: d = 3: S* (ordinate) vs M (ab- 
scissa) 



Our main findings are summarized below: 

• 4 < d < 10. The behaviour is similar to the d = 4 case which is illustrated in Fig. [l|. To 
the right of their respective peaks, the temperature, mass and entropy exhibit an infinite 
series of damped oscillations: configurations in this range are unstable to one or more 
radial modes, although the precise onset of instability depends on the ensemble In the 
canonical ensemble, at fixed red-shifted temperature, all configurations to the right of T^ax 
are unstable. In the microcanonical ensemble, at fixed mass, all configurations to the right 
of Mmax are unstable. The maximum temperature configuration always occurs at a lower 
central density than the maximum mass and entropy configuration (note that for d = 10 
they lie very close) . The numerical values given in Table |l] for d = A agree with . 

• d > 11. In this case the red-shifted temperature, mass and entropy increase monotoni- 
cally with the central density, asymptoting to their maximum values as the central density 
goes to infinity, see Fig. |2|. Since the specific heat is always positive, one would expect 
these configurations to be stable for all values of pc- To the author's knowledge, this case 
has not been studied before. Note that this critical dimension was found analytically and 
not numerically: its exact (non-integer) value is in fact just under 11, see Appendix 
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• d = 3. In this special case an exact solution exists, see Appendix We found the 
following expressions for the (dimensionless) red-shifted temperature, mass and entropy: 

M = AEi^ (3.3) 

S=^^. (3.4) 
1 + Pc/2 ^ ' 

These are plotted in Figs. ^ and ^. Here the mass is defined such that it vanishes for the 
massless BTZ black hole. The minimum mass M = — I/8G3 (M = — vr) corresponds to 
AdS-s in global coordinates. The maximum mass is M = I/24G3 (M = vr/S). This is the 
same as the range of masses found for black holes localized on an AdS^ brane Q.^ The 
limiting infinite central density solution is the same as the massless BTZ black hole except 
that there is a curvature singularity proportional to a (5-distribution supported at the origin 
(Appendix |A|). As in the 4 < d < 10 case, one would expect that configurations to the 
right of the peak of temperature (mass) in Fig. ^ are unstable to radial perturbations with 
canonical (microcanonical) boundary conditions. 

4. Comment: validity of perfect fluid approximation 

Under what conditions is this a good model of quantum fluctuations about the thermal 
AdS saddle point (that is, AdS periodically identified in imaginary time)? For a general 
(non-conformal, self-interacting, arbitrary spin) field theory in a static curved spacetime, 
the energy-momentum tensor (2.4) with the Stefan-Boltzmann relation ( |2.5| ) corresponds 



to the leading term in the high-temperature expansion of the one-loop effective action 
|1C, 11, The higher order terms will be small if the red-shifted temperature satisfies 
T ^ TZ^^ and T » f^eff-, where 7^ is a characteristic curvature radius of the spacetime 
and nieff is the effective mass of the field [11 The constant a is proportional to the 



effective number of massless degrees of freedom g: 

a = {d-l)^-'''^v{^C{d)g, where g = ub + (l - 2'^''-^^) np , (4.1) 

with ub being the number of boson spin states and np the number of fermion spin states. 
This may be obtained by integrating the d-dimensional Planck distribution over all fre- 



quencies, taking into account the degeneracy of states factor, see, for example, [13|. Note 



that we are also assuming that the characteristic curvature radius TZ is much greater than 



*Fig. ^ should be compared to Fig. 1 in although the dependence of the entropy and temperature on 
the other parameters (G3, I, a) is different for the two cases. The expressions for the mass, temperature 
and entropy in also depend on the parameter A, which is related to the positions of the two branes. 

^In general, the effective mass will contain a term involving the Ricci scalar, as well as terms involving 
the interaction potential. However, since we are expanding about the trivial ^ = background, there will 
not be any interaction terms. 
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the Planck length Ipianck so that vacuum polarization terms and higher order quantum 
gravity effects are negligible. 

In the context of the configurations studied in this paper, we can be more precise. 
If the field is massless (or if the red-shifted temperature is well above the rest mass), the 
condition for the higher order terms in the high-temperature expansion to be small becomes 

T » max{(K/)c)^/^ , r^}, (4.2) 

since for npc ^ the origin will be the place of maximum curvature. Now dip/dr > 
and condition ( p.3| ) imply 1^(0) < 0. Prom equation ( p.9D and V{0) = it follows that 
21oc(0) > ^- Equation ( |4.2| ) and the Stefan-Boltzmann relation ( |2.5D may then be used to 
show that the higher order terms are small if the dimensionless temperature satisfies 

or if the dimensionless density falls in the range 

7"' « Pc < 7'/'' , (4.4) 

where 7 = a~^{l/lpianckY~'^ ■ For / » Ipianck, as we have assumed, and if the number 
of degrees of freedom is not too large, 7 is a large number and this model is applicable 
in a wide domain.^ In that case, in 3 < d < 10, the majority of the configurations 
corresponding to the curves in Pig. || and Pig. |^ would belong to this regime, including the 
maximum temperature configuration. On the other hand, in d > 11 (Pig. §), the maximum 
temperature configuration corresponds to infinite central density and therefore higher order 
terms would come into play well before this point is reached. It would be interesting to 
understand better the significance of this. 

The same factor 7 appears if we compare the maximum temperature of thermal radia- 
tion to the Hawking-Page temperature Thp = (d— 2)/27r/: we find Tmax/Tup ~ 7^'^'^ ^ 1 if 
7^1. This is consistent, since thermal radiation heated to a temperature T > T^ax must 
collapse to a black hole with a lower free energy. Similarly, if we follow |l| and consider 
the microcanonical ensemble at fixed energy E, it may be estimated that an equilibrium 
of a black hole with thermal radiation is more probable than thermal radiation alone if 
E > El, where 

Por consistency we require Mmax > El and indeed we find Mmax/Ei ~ 7('^-3)/(2'i-3) ^ 
1, assuming d > SJ Purthermore, we may compare the maximum entropy of thermal 
radiation to the entropy of a black hole of the same mass, SpH = 27rS/t~^r^^, where 
rn is the horizon radius: we find Smax/SpH ~ 7^^/^ ^ 1, which is consistent with the 
Bekenstein bound. Note that these comparisons should not be taken too literally, since we 
should really be comparing the radiation values to the black hole values calculated at the 
one-loop level. 



The upper limits in (4.3) and (4.4) are well below the Planck scale in that case, 
''in d = 3, there exists a stable equilibrium of a BTZ black hole and thermal radiation for any E > Q. 
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5. Conclusions 



We have studied spherically symmetric equilibria of perfect fluid radiation with asymptot- 
ically AdS boundary conditions in various dimensions. As in the d = 4 case, there is a 
maximum red-shifted temperature, maximum mass and maximum entropy in all dimen- 
sions d > 3. We computed these values for 3 < d < 13. In d = 3 the analytical results 
serve as a useful toy model, although it must be kept in mind that 2 + 1 gravity has some 
special properties that do not necessarily carry over into higher dimensions. In 4 < d < 10 
the behaviour is similar to the d = 4 case studied in In d > 11 there is a qualitative 
difference in behaviour: the temperature, mass and entropy are monotonic functions of the 
central density, approaching their maximum values as the central density tends to infinity, 
with no oscillations. Whilst we are not certain of the significance of the latter result, what 
we can say is that it applies in a more general setting than that of a radiation fluid. The 
oscillatory profile of thermodynamic variables at high central densities, associated with 
modes of instability, is a feature of a wide class of stellar models: it appears in general 
relativity whenever the high density core of a star is described by a gamma-law equation of 
state [ 14 , 0, |l^ , or indeed certain polytropic equations of state [^] , as well as in Newto- 
nian theory (the isothermal sphere ||l^), and is independent of the properties of the outer 
regions ||l^. As shown in the dynamical systems analysis of [^], for an asymptotically 
gamma-law equation of state, it can be traced to the nature of the self-similar solution 
(2.14), which is a fixed point of the Tolman-Oppenheimer-Volkoff equations. In d = 4 it is 
a stable focus, resulting in oscillatory behaviour in the high density limit. Their analysis 
may be extended to d dimensions, see Appendix ^ We found that in 4 < d < 10 it is a 
stable focus, whereas in d > 11 it is a stable node, resulting in monotonic profiles, as seen 
in this paper. 

One may also look at these solutions from another perspective, which comes from con- 
sidering the connection between dynamics and thermodynamics 0]. The entropy of a given 
state is expected to measure the logarithm of the fraction of time the system spends in that 
state throughout its dynamical history ||2^. For matter configurations, this has the impli- 
cation that the stable equilibrium configurations should correspond to local maxima of the 
entropy at fixed energy [Q]. The relevance of this is as follows. Firstly, it may be shown 
that all static spherically symmetric asymptotically AdS configurations of a radiation fluid, 
i.e. the solutions to the Tolman-Oppenheimer-Volkoff equations ( |2.6D are extrema of the 
total entropy ( |2.12| ).^ The converse statement, that the solutions to equations ( |2.6D are 
the only extrema of the entropy, would hold only if it can be assumed that all extrema 
are spherically symmetric and that all fluid spacetimes of the type considered contain a 
maximal hypersurface @. Existence of maximal hypersurfaces in asymptotically flat |21| 
and asymptotically AdS |22| spacetimes has been proved subject to certain conditions. The 
long-standing conjecture that a static perfect fluid star in d = 4 is necessarily spherically 
symmetric appears to have been proved in |^] for physically reasonable equations of state 
in asymptotically flat space. It is not known if a result of this kind exists in asymptoti- 



*This was shown in Q for the A = 0, d = 4 case of radiation in a box of fixed radius, but it extends to 
the asymptotically AdS case. 



-10- 



cally AdS space, in arbitrary dimension. Note that in the asymptotically AdS case one 
can consider more general conformal structures at infinity, the so-called asymptotically 
locally AdS spacetimes [^]. Thus in this paper we only dealt with a restricted class of 
solutions. Secondly, [Q] showed that in the A = 0, d = 4 case of radiation in a box, the 
criterion for thermodynamic stability - that the extremum is a strict local maximum of 
the entropy - is equivalent to the condition for dynamical stability under linear radial per- 
turbations. This stability argument would be expected to generalize to the asymptotically 
AdS case in d dimensions, although it would be useful to carry out a detailed perturbation 
analysis. A similar formulation should also be possible in the canonical ensemble, where 
stable equilibria should correspond to local minima of the free energy at fixed red-shifted 
temperature. 

Finally, it would be interesting to consider these solutions in the context of the Karch- 
Randall model, where one has an asymptotically AdSd-i brane in an AdS^ bulk |25, 26|. 
According to AdS /C FT, classical dynamics in the bulk is dual to the quantum dynamics 
of a CFT living on the brane, in the planar limit, coupled to {d — l)-dimensional Einstein 
gravity. The strongly-coupled CFT has a UV cutoff ^ 1/L and an IR cutoff ~ 1//, where 
L and I are respectively the bulk and brane AdS lengths. This theory has been interpreted 
as a defect CFT although there are a number of issues still to be resolved. In the 
original Karch-Randall model, the bulk volume is infinite and the zero mode graviton is 
not normalizable, nevertheless for I ^ L there is an ultralight graviton with a tiny mass 
m? ^ L'^^^/l'^^^ and (d— l)-dimensional gravity is effectively reproduced at distances r > L 
on the brane, although its mass, and the extra dimension, would become apparent at very 
large distances r > l'^~'^/L'^~^ [^1]-^ However, if one introduces into the bulk a second 
positive tension AdSd-i brane [^, one recovers the zero mode and (d— 1) dimensional 
gravity is valid all the way to infinity. One now has two gravitons: the massless one and 
the ultralight one, although the ultralight one decouples in the limit that the second brane 
approaches the turn- around point of the warp factor |29|. Now if the metric on one of 
the branes is given exactly by one of the thermal radiation configurations studied in this 
paper, one may wonder what happens in the bulk? In this context, it is interesting to note 
that in the d = A version of this model, with two AdS^, branes in an AdS4, bulk, the bulk 
metric is the AdS C-metric studied in Q and the other brane contains a localized black 
hole (Appendix ^). It would be interesting to see if this property is just a special feature 
of this lower dimensional model or whether there is a regime where it would hold in higher 
dimensions. We note that the issue of the Hawking-Page transition in this setup has been 
addressed in |^^, although only in the limit //L — > oo where matter effects can be ignored. 



Note 

After this work first appeared, two other papers came out, one by John Hammersley 
and one by Pierre-Henri Chavanis |41], who have been independently working on similar 
problems. Specifically, [^] studies self-gravitating radiation in d-dimensional AdS and also 

^For L < r < I the {d — l)-dimensional graviton is a composite made out of the ultrahght mode and 
heavier Kaluza-KIein modes, whereas for I ^ r < / L'^^'^ it is just the ultralight mode [M. 
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finds the change in behaviour from oscillatory in 4 < < 10 to monotonic in d > 11 (note 
that the d = 5 case had previously appeared in in a study of bulk-cone singularities in 
AdS/CFT). He points out the possibility that (besides a possible link to string theory!) 
this could be related to a similar dimensional transition found in a cosmological setting 
1 37, where the analysis of Belinsky, Khalatnikov and Lifshitz (BKL) ||3^ was extended 
to arbitrary dimension, with the finding that the chaotic oscillatory dynamics of the grav- 
itational field close to a spacelike singularity for 4 < d < 10 is replaced by a monotonic 
power law time dependence for d > 11. It would be interesting to explore this connection 
further. 

The work of |41| is a comprehensive study of the behaviour of general relativistic perfect 
fluid stars with a linear equation of state in d dimensional asymptotically flat space (in a 
box), including the radiation case, extending the analysis of |17| to arbitrary dimension 
(the Newtonian case had previously been treated in |35|). He also finds that for d > 11 
the oscillatory behaviour disappears. In fact the results of Sec. 6.4 in |41] are consistent 
with the findings noted in Appendix ^ of this paper and provide an alternative derivation 
of the asymptotic behaviour of A = linear perfect fluid solutions in different dimensions. 
The paper |41] also includes a detailed stability analysis of the different regimes, with 
the suggestion that for these solutions thermodynamic stability is equivalent to nonlinear 
dynamical stability. The d = 3 case, treated in Appendix^ of this paper, is also analyzed, 
although the author suggests that these solutions have finite radial extent whereas, as noted 
in Appendix the finite radial coordinate at which the density vanishes is at infinite proper 
distance from any other point in the spacetime. 
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A. Exact solution in d = 3 

In this appendix we derive exact solutions for regular equilibrium configurations of thermal 
radiation in d = 3, both for A < 0, and for A = 0. Garcia and Campuzano have 
obtained all static circularly symmetric perfect fluid solutions in 2 + 1 dimensions, and 
the solutions below belong to one of the classes presented in that paper. We found it 
interesting, however, that the case of a radiation fluid with A < may also be obtained 
starting from the AdS C-metric studied by Emparan, Horowitz and Myers [^, describing 
black holes accelerating in AdS^^. This metric, which is a static axially symmetric vacuum 
solution to the four-dimensional Einstein equations with negative cosmological constant, 
takes the form 

ds- = I 

A^{x — y)"^ 



H{y)de 



+ 



dx' 



Hiy) Gix 



+ G{x)d(t)'^ 



(A.l) 
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where 



H{y) = -A + ky" - 2// V 

G{x) = l + kx^ - 2fiAx^ , (A.2) 

with A>0, y4>0, /x>0 and k = —1, 0, +1. Points with x = y correspond to the 
boundary of the asymptotically AdS4^ geometry, therefore the solution is restricted to the 
region y < x. The coordinate <j) is an angle and each zero of G{x) corresponds to an axis 
for the rotation symmetry, where we need G{x) > to preserve the Lorentzian signature. 
For /i > 0, G{x) has only one positive root X2 and in order to avoid a conical singularity 
at this point, the periodicity of the angular coordinate cp is set to be A(j) = An/IG' {x2)\- 
The smallest zero of H{y) defines the black hole horizon. See for further details. The 
authors of considered putting two positive tension asymptotically AdS^ branes in this 
spacetime: one at x = 0, the black hole brane, and a second non-singular brane at y = 0. 
This is a lower dimensional version of the two-brane Karch-Randall model where, at least 
in a certain range of parameter space, bulk gravity in AdS/i is expected to reproduce AdS^ 
gravity on each brane, plus quantum corrections coming from a CFT. The different values 
of k correspond to different slicings of AdS^.^^ Now, if there is a quantum corrected black 
hole on the first brane, then the other non-singular brane at y = should correspond 
to thermal radiation in AdS^ at constant red-shifted temperature. This may verified by 
examining the induced metric on the surface y = 0, 

dx'^ 



A^ 



x^ 



G{x) 



(A.3) 



where x is restricted to lie in the range < x < X2. The Einstein tensor for this 3- metric 
is 

G'i, = A^dmg{l, 1, 1) + ^JLA^ x^d\a,g{-2, 1, 1), (A.4) 
which describes a. p = p/2 perfect fluid with density 

(A.5) 



l^A^x-' 



SnGs 

in a spacetime with cosmological constant A = —l/l"^ = —A^. To convert this to the 
standard BTZ form, that is 

ds^ = -N^dt'^ + -SG-srnj dr^ + r'^dO'^, < < 2tt , (A.6) 
where N(r) is the lapse function, change variables to the following: 



A\G'{X2)\ 

|G'(X2)L 



(A.7) 



In the /X = case when the bulk metric (A.l) is locally AdS4, k = —1 corresponds to global coordinates 
on the AdSa slices, whilst k — and k — 1 give the metrics of the massless and massive BTZ black holes 
respectively tel. 
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Observe that x = X2 corresponds to r = 0, and x = corresponds to r = oo. It is not 
possible to express the functions m, p and A*" in a simple form in terms of r, nevertheless, 
it is possible obtain explicit expressions for the red-shifted temperature, mass and entropy. 
In terms of the dimensionless variables (|3.1|), we find 



T 
M 

S 



5 

Pc 

1 + pj2 

^{Pc - 1) 

(1 + 

1 + pj2 ■ 



Note that in the metric ( [A.3| ), the different values k = —1, k = and k = 1 are mapped 
to the ranges < pc < 1, Pc = ^ and pc > 1 respectively. Consider now the limit pc ^ oo. 



Equations ( [A.2| ) and ( A. 5 ) require pA and X2 ~ {2pA) oo. The limiting solution 
has p = and m = for all r > which is the same as the massless BTZ black hole. 
However, there is a curvature singularity proportional to a (5-distribution supported at the 
origin. 

The A = case may be obtained from the metric ( |A.3| ) by making the transformation 
X = Ax, (f) = A^^cj), G = A^G and taking the limit ^ ^ 0. Dropping the hats, this gives 
the 3-metric 

1 r „ ^-7-2 „1 



G[x) 



where G{x) = kx'^ — 2px^ and the period of cj) is An /\G' {x2)\- The Einstein tensor for this 
metric is 

G(; = ^x3diag(-2,1,1), (A.9) 
which describes a p = p/2 perfect fluid with density 

2px^ 

and A = 0. Making the choice k = 1 (this is the only consistent choice), we restrict to 



< X < X2 = l/2p so that G{x) > 0). We may compare the metric (|A.8D to a more 
standard form for A = static solutions in 2 + 1 gravity: 



ds^ = -N^dr + {l-4G3my^dr^ + r^d9^, < 6^ < 27r , (A.ll) 

where N{r) is the lapse function. Now in three dimensions the Weyl tensor vanishes so the 
spacetime outside a point mass is locally flat [32|: the point mass only affects the spacetime 
globally, making it into a cone, the value of the mass being given by the conical deficit 
angle at infinity, S^o = SnG^m. Note that this differs from the mass defined by the metric 
( |A.6| ) with A = 0. In particular m = l/AG^ in the normalization of (A.ll) corresponds to 



m, = in the A = version of (|A.6|). Note also that in the case under consideration, we 



'^This differs from the usual form |32| ^] by a trivial change of variables. 
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are dealing not with point sources but with smooth perfect fluid configurations. Changing 
variables to 



\G'{X2)\ X 



\G'{X2) 



we see that x — X2 corresponds to r — 0, and x — corresponds to r — i^max — 

4// > 0. 

The radial coordinate is bounded: < r < rmax- After some algebra we obtain the 
following expressions for the mass, density and lapse function of the regular radiation fluid 
configurations with A = 0: 



1 

m 



4G3 L 
P = P> 



1 - (1 - 27rG3Pcr^ 



l-27rG3Pcr')' 



iV = iVo (1 - 2^G3Pcr') ' , (A.13) 

where pc = (327rG3/u^)~^ is the central density. The mass varies from m = at r = 
to m ^ I/4G3 as the radius r — > r^aax- This is exactly the range of masses one would 
expect in 2 + 1 gravity with zero cosmological constant with smooth matter sources p3| . 
It may be verified that the proper distance between two points in the spacetime diverges 
if the radial coordinate of one of them approaches r^ax '■ thus r = r^ax lies at infinity and 
the solutions are unbounded. There are no conical singularities: the conical deficit angle 
increases smoothly from zero at r = and tends to 27r as r ^ rmax- One may visualize 



the geometry represented by the the spatial part of the metric (A. 11) by embedding it 
into a 3-dimensional flat Euclidean space {r,9,z): one ends up with a semi-infinite cigar, 
asymptoting to a cylinder at infinity.^^ The cylinder is the d = 3 analogue of the singular 



self-similar solution ( 2.14 ): it is the geometry of a point particle at the origin of mass 
M = I/4G3 (in particular, p = for all r > 0) (3^]. As in the d > 3 case, this solution also 
corresponds to the p^^ 00 limit of the A < radiation fluid configurations for r/l ^ 0. In 
that case, the second of equations (|2.6| ) implies that p = for all r > 0. This explains our 
previous observation, that the pc ^ 00 limit of the A < radiation fluid configurations is 
the same as a massless BTZ black hole for r > 0. 



B. Phase plane analysis of A = scale-invariant system 

In this appendix we outline a dynamical systems analysis of the d-dimensional Tolman- 
Oppenheimer-Volkoff equations in the A = case, which yields an exact (non-integer) 
value for the critical dimension at which the oscillatory behaviour ceases. For generality, 
we assume a linear ("gamma-law") equation of state p = qp, with < g < 1: this leads to 

^■^As noted in |^^, these unbounded solutions are marginally stable (as in the d > 4 case). However, one 
may enclose the fluid in a box of finite proper radius and match it on the outside to a cone to obtain stable 
asymptotically conical solutions. 
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the equations being scale invariant, which greatly facilitates the analysis. We also restrict 
attention to d > 4 (the d = 3 case must be treated separately, we do not describe it here). 
In d = 4, analogous phase portraits have been considered in |Q, 17| using non-compact 
variables. We follow the formulation introduced in |19|, in which the state space of the 
dynamical system is compact, i.e. the boundaries are included. 

The Tolman-Oppenheimer-Volkoff equations in d dimensions with A = read: 



dm 
dr 
dp 
dr 



id 



{d 



- 3)Kmp f-. P 



1 + 



d-l 



P 



{d — 3)m 



2Km 



{d - 2)Sr^-3 



(B.l) 



where the metric is as before, i.e. (2.1) with A = 0. They have been studied before in 
where a higher-dimensional version of Buchdahl's limit was derived. We are interested in 
the linear and homogeneous equation of state p = qp which is the only equation of state 
that leaves (B.l) invariant under the scale transformation 



Cr, p^C-^p, p^C-^p, 



m 



d-3 



m . 



(B.2) 



where C G 



H]. In that case, introducing the dimensionless variables 



13 



d-l 



P 



Km 



m 



1 



2 Km 



p{d- 2)i;r°'-3 \^ (^d- 2)Sr^-3 



(B.3) 



equations ( B.l ) are transformed (cf. |19|) to the following autonomous system: 

du 

= -u{l-d + u + {d-3 + qu){l + q)v) 

d^ 
dv 

— = —V {d — 3 — u + 2q {d — 3 — u) v) , 
dt, 



(B.4) 



where ^ = Inr. Note that if m and p are positive and the spacetime is static, u and v lie in 



the ranges u > 0, v > 0. In that case, following |19], we change from {u,v) to the bounded 
variables {U, V) G (0, 1)^ 

U = -^, V = -^, (B.5) 



l + u 

and introduce a new independent variable A, 

d\ _ 
d^~ 



1 + v 



1 



{l-U){l-V) 



(B.6) 



^^These may be regarded as a relativistic analogue of the Milne variables, homology invariant variables in 
the Newtonian theory. If the speed of light c is put back in, we have for a linear equation of state p = qpc^ , 
where p is now the rest frame mass density. The Newtonian limit is given by c — > oo, g ^ with qc^ 
constant. The Einstein equations then reduce to an Emden equation describing an isothermal sphere in 
Newtonian gravity . Note that in the Newtonian case a polytropic equation of state also preserves scale 
invariance 
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2d 



whence equations ( [B.4D yield the following dynamical system: 
dU_ 

dV 

As discussed in [jl9|, the asymptotic properties of a dynamical system are determined by 
the boundaries, therefore it is essential to consider the compactified state space [0,1] 2. This 
requires the system to be on [0, 1]^, which is clearly the case for the polynomial system 



U{l-U){d-l-dU 
V{l-V)il + i2q-l)V) (3 



- 4 + {d - 3) q)V + {2d - 3 + {d 
d+{d-2)U). 



4) q) UV) 

(B.7) 



(B.7). 



The idea is to solve the dynamical system (B/7) for U{X), V{X) on the state space 
[0,1]^ and then to obtain the physical perfect fluid configurations m(r), p{r) one inverts 
the transformations ( p.3| ), ( |B.5| ), see ]T9| . 

A pivotal role in any dynamical systems analysis is played by the fixed points. The 
fixed points of ( |B.7D in the unit square [0, 1]^ and their associated eigenvalues are listed in 
Table |2|, where the parameters ui and uj2 are given by 
{d-3){l+q) {d-3 + {d-l)q) 



002 



2{d-2) {d-l + {d-3)q{2 + q)) 

{d - 3) (1 + q) + q)^d'^ - 2(1 + q){7 + 5q)d + 25^2 j_ 22q + 33 





2{d 


- 2) (d - 1 + (d - 


-3)qi2 + q)) 


Fixed point 


U 


V 


Eigenvalues 


Ti 








[d-l) , -{d-3) 




d-l 
d 





(d-l) 2 

d ' d 


T3 


d-3 
d-2 


2 

d-l+(d-3)g{2+g) 


—iOi ± UJ2 




1 





1 , 1 


n 


1 


1 


q{l+q) , -2q 







1 


-{d-3){l + q) , 2{d-3)q 



(B. 



Table 2: Fixed points of (B.7) in the unit square [0, 1]^ and associated eigenvalues, assuming 
< g < 1 and d > 4. 

The global dynamics of the system is determined by its past and future limit sets. 



For the system (B/7) it takes a particularly simple form: it may be shown (cf. |]1^]) that 
orbits in the interior of the unit cube originate either from the fixed point T2 or from the 
fixed point and converge to the fixed point T-^}^ The central role is played by the 

^"^The Hartman-Grobman theorem states that the eigenvalues and eigenvectors of the hnearization of 
a dynamical system about a hyperbolic fixed point completely characterize the flow of the system in the 
neighbourhood of this point, that is, the flow of the full non-linear system is topologically equivalent to the 
flow of its linearization in this neighbourhood. 

^^This follows from local analysis and the monotonicity principle, whereby the existence of a monotone 
function on the state space precludes periodic or recurrent orbits, see |l9j. It may be verified that a 
monotone function is given by Z = {qhi + ki)'^{uv)^^''-{l + 2qvY''-~^ where /ii = (1 + q){d — 3 + qu)v, 
ki = d - 3 + {d - l)q and pi = V/ ((1 + q){d-S + {d+ l)q), cf. ll| for d = 4. 
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fixed point T2, corresponding to Minkowski space in splierically symmetric form, and by 
the fixed point T3, wliich represents tlie singular self-similar solution (cf. eqn. ( p.l4D ) 

p = ar'^, m = aT.r'^^'^ (d - 3 where a = — '-^7: -. (B.9 

Scale invariance implies that the entire set of solutions with a regular centre is repre- 
sented |19| in the state space by a single orbit originating from T2 (as A —00 or r — > 0) 
and converging to T3 as (as A — > 00 or r — > 00).^^ The exact form of this orbit must 
be found numerically, however, one may write down asymptotic approximations for it by 
linearizing the dynamical system about the points T2 and T3: these are determined by the 
corresponding eigenvalues (and associated eigenvectors) in Table ||. It follows from ( |B.g| ) 
that the regular perfect fluid configurations with A = have infinite mass and are thus 
not asymptotically flat: to obtain finite mass perfect fluid solutions one must confine the 
radiation to an unphysical box. Scale invariance then implies that the regular orbit admits 
a dual interpretation Q, whereby different initial segments (equivalently end points) of 
the orbit are interpreted as representing distinct regular equilibria of a perfect fluid in a 
box of flxed radius R, that may be parametrized by their central density pc (a monotonic 
function of A), so that T2 represents the solution with pc = 0, i.e. an empty box, and T3 is 
the limiting solution as pc ^ 00. 

The crucial fact is that the behaviour of the regular orbit as A — s- cxd, i.e. in the 
neighbourhood of the fixed point T3, reflects itself in the high-density behaviour of various 
asymptotic quantities, such as the mass, the red-shifted temperature or the entropy, of a 
linear (i.e. p = qp) perfect fluid in a box. Indeed, it was shown in |jl9| that the equation 
of state may be quite different away from the centre, as long as it is still asymptotically 
linear at high densities. They considered a fairly general class of equations of state that 
are asymptotically linear at high densities and asymptotically polytropic at low densities 
and proved that the mass-radius diagram of these stellar models has a spiral structure (for 
low enough polytropic index the radius is flnite without the need for a box).^'' The case 
treated in this paper, where the linear equation of state was modifled by adding a negative 
cosmological constant, giving finite mass unbounded configurations without the need for a 
box, may be treated similarly. 

The behaviour of the regular orbit near T3 is determined by the eigenvalues — ±^2; 
where uJi and uj2 are given by equation ( |B.8| ). These depend on the value of q and on 
the dimension d. For < (7 < 1 and 4 < d < 10, is always imaginary, and the fixed 
point r3 is a stable focus (since uJi is positive). This leads to oscillatory behaviour, e.g. 
the spiral structure theorem in |19| or the 4 < (i < 10 case in this paper. However, for 
a given < g < 1, there is a critical dimension dcrit above which u)2 is real and strictly 



^^The one-parameter set of orbits originating from the fixed point Ti may be continued to the negative 
mass domain it < 0, n < and describes solutions with a negative mass singularity surrounded by positive 
energy density that enter the state space when they have gained enough mass fl^ . 

^'^The presence of another, non-scale-invariant, variable in the equations means one has to work with a 
3-dimensional dynamical system, i.e. the state space is enlarged to a cube |l^]. The key is again scale 
invariance, which is now only present asymptotically |19|]. 

^*For the case of positive cosmological constant, see jl8|. 
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less than uJi, i.e. T3 is a stable node and monotonic behaviour emerges instead. It may 
be verified that the critical dimension is always in the range 10 < d < 11. The case of 
perfect fluid radiation {q = l/(d — 1)) treated in this paper gives a critical dimension 
dcrit = 10.96404372... Curiously, a stiff fluid {q = 1) gives dcrit = 10, and pressureless 
dust {q = 0) gives dcrit = 11. The last result also corresponds to the Newtonian limit (the 
isothermal sphere) and has been obtained before in [ p5[ . 
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